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SUMMUtT 


V 

The  problem  of  wave  propagation  In  a  medium  with  fluctuating  refractive 
Index  la  considered,  this  problem  being  of  interest  fro*  the  standpoint  of 
the  theory  of  twinkling  of  stars  and  atmospheric  acoustics.  The  problem  is 
reduced  to  linearised  equations  satisfied  by  the  phase  and  lofpxlthsdc 
amplitude  of  the  propagating  ware.  Limiting  cases  are  investigated  and 
particular  exasq^les  are  considered  to  indicate  the  listlts  of  validity  of 
the  jreometric  optics  treatment  and  the  nature  of  the  corrections  required 
to  take  account  of  diffraction  affects.  ,  _ 
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Small  atmospheric  inhcmofenieti.s  ari.lng  fro.  turtulenc.  vurt  a 
profound  lnflu.no.  oo  propagation  cf  scxtnd  and  light,  causing  the  characfr- 
l.tlc  fluctuations  of  ti  i  associate  fi.ld.  t.  k.  Krasilnikov  has 

investigated  experimentally  phase  and  amplitude  fluctuations  of  sound  wavea 
propagating  through  the  atmosphere.  Similar  fluctuations  of  light  waver  ara 
exhibited,  for  instance,  by  the  phenomenon  of  twinkling  of  stare,  which  are 
evidently  cloeely  connected  with  turtxilent  fluctuations  of  atmospheric 
temperature.^’  **3* 

Reference*  1,  3  and  8  contain  some  computations  of  the  phase  and 
amplitude  fluctuations  of  a  wave  propagating  in  a  turbulent  medium  for 
oblique  incidence.  However,  this  work  is  based  on  the  geometric  optica 
( acoustics)  approximation,  which  evidently  in  certain  cases  is  inadequate 
and  gives  rise  to  a  discrepancy  between  computed  and  experimental  values. 

For  instance  in  the  geometric  optics  approximation  the  amplitude  fluctu¬ 
ations  are  proportional  to  the  3/2  Fewer  of  the  (turbulent)  layer  thickness, 
wnereas  observations  generally  indicate  a  considerably  lower  order  of  grewth. 

In  the  present  paper  an  attempt  is  made  to  consider  the  jrcbleir.  of 
amplitude  and  phase  fluctuations  cf  a  scalar  wave  field  on  a  more  general 
basis  and  by  means  of  more  accurate  approximations  taking  into  account 
diffraction  effects,  to  thus  clarify  the  tints  <_r  validity  cf  the  geometric 

optics  approximation. 

1.  Derivation  of  the  Baalc  Squat! one. 

Let  us  consider  a  scalar  wave  field  riven  by  V(x,  y»  a,  t)  ■  >  (r,t) 
which  satisfies  th«»  equation 


d  4  = 

f  H*  ' 


(1) 


*A  theoretical  armlysls  of  turbulent  fluctuations  cf  t:;e  atmospheric 
temperature  and  certain  related  experimental  data  can  be  foun)  in  references 


L 5  'll- 
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where  the  propagation-velocity  c  of  the  wave  is  a  function  of  position, 
the  time  dependence  of  c  will  be  neglected,  assuming  that  the  characteristics 
of  the  medium  are  slowly- varying  in  time  (i.e.,  with  respect  to  the  carrier 
frequency).  Further  it  will  be  supt oaed  that  the  deviations  of  o  from  some 
nean  value  e  : 


c 


(2) 


are  sswll  everywhere  in  the  modium  and  for  all  time.  The  refractive  index 
n  fluctuates  around  unity 


n  ■  1  ♦  p,  (n  <<  1) 


U') 


(in  this  case  it  Is  possible  to  speak  of  a  slightly  inhomogeneous  Md  us), 

Let  us  suppose  that  the  half-space  x  <  0  is  homogeneous  (i.e.,  n  ■  1 
for  x  <  0)  and  that  plane  waves  are  incident  at  the  plane  x  ■  0  from  the 
direction  jf  the  negative  x-axis.  The  problem  consists  of  determining  the 
cliaracterietlcs  of  the  wave  field  in  the  plane  of  observation  x  •  L.  It  1« 
clear  that  this  problem  is  closely  related  to  the  problem  of  diffraction  of 
light  by  ultrasonic  waves,  which  has  been  Investigated  by  3.  M.  Rytovt?! 
but  however  differs  in  that  we  canr.ot  assume  periodic  variation  of  tue 
refractive  index. 

Coneidur  the  case  when  <p  r,  t)  is  a  mt ncchroewtl c  wave  of  frequency  wj 


■f  (r,  t)  -  A(r )  exp  <  1 


<**  -  S(r) 


(3) 


For  the  optical  problem  (e.g.,  twinkling  of  stare)  the  scalar  aquations 


apply  if  polarisation  effects  may  be  neglected. 
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vlth  amplitude  A(r)  and  phase  S(r).  It  is  appropriate  to  introduce  a 
coiflex  function 

y  (r)  -  3  ♦  1  In  -p-  , 


U) 


ao  that 


<p(r,  t)  -  A  exp 


«t  -  +(r) 


(5) 


where  Aq  1b  the  amplitude  of  the  (unperturbed)  Incident  wave.  We  have  thus 
passed  from  equation  (1)  to  the  elkonal  equation  (s_lc) 

.  i  i 


(v+)  +  i dt  =  L ”, 


t,-.  ^ 


(6) 


In  the  case  of  a  homogeneous  medium  we  would  have  In  place  of  equation 
(6)  the  equation 


(6  ) 


(Y"  is  the  equivalent  of  V'  for  a  homogeneous  medium).  Subtracting  (6/)  from 
(6)  yields 

2(v+-vt'>  +  ■  at'  =  _><*  +  l/ii-  (?r>‘] 


where 


^  --t-'K, 


py/  _  _  py- 


f 

Having  asmimed  that  la  of  the  order  of  p,  ore  has  a  basis  for 

2 

neglecting  terms  of  order  p  In  tbs  square  brackets  cn  the  right.  We  then 
obtain  a  linear  equation  in 
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2  V^'-Vt0  -h  i  A  +  '  =  */, 


(7) 


which  will  serve  as  a  baala  for  the  further  analysis.  Similar  linearised 
elkon&l  equations  were  employed  for  the  first  time  in  the  wcidc  of  3.  M. 
Rytovk?^  in  1937. 

/ 

We  note  that  ■  S  -  ia  the  fluctuating  (perturbeu)  phase  of 

the  ware  brought  about  by  fluctuations  of  the  refractive  index  n,  and 
Y  ■  In  (*A  )  logarithm  of  the  ratio  of  the  amplitudes  of 

perturbed  ami  incident  fields.  Because  in  the  derivation  of  equation  (7) 

/ 

we  have  supposed  only  that  V7 /r.  <  ■»  1,  it,  is  clear  that  no  requirement 
of  smallness  of  phase  and  asiplltude  perturbations  (with  respect  to  n  and 
lo,  respectively)  has  been  imposed.  What  is  required  is  that  the  chance 
of  phase  and  logarithmic  Mqjlltude  over  one  wavelength  (in  an  arbitrary 
direction)  be  eufficiently  ssudl.*  This  condition  is  satisfied  for  Instance 
by  the  refraction  of  the  primary  ray  at  small  angles  or  by  diffuse  radiation 
provided  the  energy  of  the  diffuse  component  is  email  in  cosparison  with 
the  energy  f  the  primary  field. 

If  one  substitutes  in  (7)  the  function  ^  •  k  x,  c or r 6 spending  to  a 

o  0 

plane  wave  pr  pagating  in  tne  direction  cf  the  positive  x-axla,  there 
results  the  equation 


,4V 


or 


(8) 


Fr  t,  this  equation  it  is  possihle  to  obtain  the  geoewtrlc  optics 


We  note  that  th<*  classical  method  -f  solving  the  wave  equation  (1)  by 
means  of  f.  methoi  of  f  <*rturbati  on  t:.eor\,  which  is  often  applied  in 
scatterinr  problems,  becomes  lnapjlicafcle  at  large  I.  because  <f  "phase 
pillnr  up"  to  values  ccmj^rable  with  v  (cf.  [lo]). 


approximation  in  the  way  used  by  V.  A.  Krasilnikov**..  Let  us  introduce 

the  optical  path  lenfth  %  -  S/A#  and  *•  -ate  real  and  imaginary  parts  of 

? 

(8),  having  previously  divided  the  whole  equation  through  by  2k0  . 


(9) 


JA.  (A /Ae) 


=  O. 


(10) 


Assuming  that  for  kQ  — p  oe  the  quantities  ®  and  In  (A/A^)  tend 
toward  finite  lialts,  cna  obtains  in  the  lisiit  the  equations  of  V.  A. 
Krasilnikov! 

j?'r  ,, 

that  is 


t.'  -  r  j 
b  -  j  h  d* 


(11 


j  A&  Mi 

A  -v 


(12) 


or,  substituting  (11)  into  (12), 

-  j  i  ...(uj-i  (a -.)(p 


[<"  £  ] 

U  -i 


'i  :  “(c) -'p 

L 


~  J  tL-^AfU  dir 

n  ^  / 


(13) 


^  ^  *■ 

where  A.  -  —  ■*-  is  the  "transverse*  Laplaoe  operator.  Bata  in  lag  tbs 

-i  jy* 

diffraction  Urn  —  [jU  A/A0)  i“  equation  (9)  aakes  possible  aera 

general  results  by  approxiaatiog  tbs  wave  effects  and  thus  peraita  one  to 
obtain  estiaates  of  the  Halts  of  Krasilnikov's  theory. 

A  solution  to  equation  (6)  will  be  sought  of  tbs  fora 

y  =  ey( iko*)  ”■  (U) 

Substitution  of  (H)  into  (8)  shows  that  w  satlsfiss  tbs  InboaDgeteous 
wave  equation 

A  kv  +  *o  ^  =  “  f(*jy.‘x )  (15) 

where 

j  -  2  ,/<  4,  e»j>(-«‘ >.«), 

It  is  well  known  that  the  solution  to  equation  (15)  oan  be  written  in 

s 

the  fora 

«<«  *  rn  III  ’  -^2^  f  ("UK, ,  d‘) 

where  T  is  the  dona  in  of  integration  and 

r  *  J/ ( >  -  t  )2 ' y  7  f-i~('r-C) 

is  the  dlstanoe  between  the  point  of  observation  P(x,y,i)  and  the  variable 
point  M(  <*  If  as  now  pass  back  to  the  function  S'  ,  we  find 


e 

The  question  of  what  conditions  to  lapose  on  w  to  exclude  singular 
solutions  will  not  detain  us  here.  In  conarete  examples  it  is  not  difficult 
to  verify  that  other  solutions  (which  oontaln  terns  of  *advanoed  potential" 
type)  lead  to  physically  absurd  expressions  for  ^ 
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x  T^‘  !  f  ("eht'lr  (—*>]} 


Thus  we  are  successful  in  obtaining  an  exact  solution  to  equation  (8), 
the  course  function  being  of  the  fora 

/  e*p5-‘tolr-/j,-f}W 


For  the  further  investigation  we  restrict  ourselves  to  the  cass  where 
the  wavelength,  though  finite,  is  several  ties a  aealler  than  the  oharaoter- 
lstlo  dlnenslons  of  the  perturbation,  that  is,  to  Fresnel  diffraction.  In 
this  oass  one  naturally  assuas a  that  the  field  at  the  point  (L,  0,  0)  results 
only  free  tbs  contributions  of  the  esdiua  inside  a  narrow  oons  with  vertex 
angle  l  <  <  l  and  vertex  at  the  eouroe |  in  other  words  it  la  possible  to 
neglect  the  scattering  et  large  angles  (in  particular,  the  backward  scattering). 


Setting 


■*  = 


f(r  +  U  -  O' 


one  sees  th*  our  assunptlons  aean  that  the  essential  contribution  to  the 


integral  (17)  coaes  fro-  that  region  for  whioh 


*-  S'*  o 


i.s.,  p/%  ~  £  <4  1 


It  then  follows  that 


(>-t)  + 


2  (*  -  \) 


* 

Adaisslbllity  of  values  of  (  can  be  decided  by  scans  of  the  scattering 
lodicatrlx.  For  the  oaas  L/X  >  1  and  weak  perturbations  this  can  be  ooaputed 
by  ths  act  bod  of  K.  S.  Shifrln. 


1 
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If  one  replaces  l/r  in  tbs  Integrand  of  (17)  by  and  r-(x-  <*  )  by 

the  aaoond  order  terms  of  aquation  (18),  than  in  the  raflon  of  important 
oontributiona  to  tbs  flald  at  (L,  0,  0)  ona  has  approximately 


Je'f>  i~,Jroir-  etp[  aCn-f)  J 

i 

Tha  expression  for  ^  (x,y, s)  can  than  ba  written  approximately  a a 


(19) 


vtw)  -  £  If  / 

T' 


(20) 


where  T'  la  tha  raflon  such  that  x-  <£  >  0  (l.a.,  altuatad  between  tha  aouroa 
of  tha  wavs  and  tba  point  of  observation).  We  raaark  that  formula  (20) 

lonatltutaa  tbs  axaot  solution  to  tba  aquation 

^yr-'  .  /  '  d 


2  k 


+  I 


t  d’'v  !  -  ?  L 

(jp  ~  ^k". 


which  differs  froa  (8)  only  in  that  tba  *transve-ee*  Laplaoe  operator 

J 

appears  instead  of  A'V  (l.e.,  tba  term  in  d  la  omitted ) .  The  offset 

0  * 


of  neglecting  tba  tara  ^  has  an  offset  on  tha  solution  similar  to  that 

3  2 

for  tha  purely  periodic  perturbation  Justified  by  S.  M.  J^ytov  j  9  j  • 

”W- 

is  a  means  of  coupe  ring  tha  results  given  by  tba  formula  ( 20 )  with  those 
following  from  tha  geometric  optics  approximation  wa  compote  tha  phase  aad 
aaplitude  perturbations  at  tha  point  (L,  0,  0)  for  tbs  following  lncraaant 
to  the  refractive  lndaxt 


i  \ 

“  =  '■/  (~  ~TF~W‘y 


(21) 


We  can  now  extend  formally  the  Integration  over  tha  whole  spaoa  T* ,  not 
only  in  tha  interior  of  tha  oontrol  cone,  because  outside  this  oone  tha 
rouroe  funotlon  (both  approximate  and  praclae;  rapidly  oscillates,  so  that 
integration  over  the  space  outside  the  cone  yields  negligible  coctributior* 
for  sufficiently  smooth  variation  of  p(x,y,»). 
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» 


(The  ecale  t  evidently  characterises  the  transverse  dimension  of  the  perturba¬ 
tion.  ) 

Performing  the  integration  in  f|  and  in  (20)  yields 

•rfL  o  o\  =  fJz+lk-ti/t* 

JoTTu^T77TTr/,M  L 

In  view  of  (21)  thla  expression  oan  be  put  in  i  forms 


L  D  c)  — 
)  /  ' 


k 


s _ 


i  t  (L-  |jV(*;  /*) 


'  f 


The  phaee  fluctuatlone  are  given  by  the  first  summand  on  the  right-hand 
eide  of  this  expression  and  the  fluctuations  of  the  logarithmic  amplitude  by 
the  second  summand.  This  result  differs  from  formulae  (11)  and  (13)  obtained 
under  the  geometric  optics  approximation  by  the  presence  in  the  Integrand  of 
tbs  factor  X/  Q  1  ■»  (L  -  £  )2/k2t*)J  .  It  follows  that  ths  geometric  option 
formulae  give  the  correct  result  for  the  case  L/(kQi2)<<  1.  Introducing  the 
wavelength  X  <=  2*/ k  enables  one  to  write  the  criterion  Just  given  in  the  form 


(22) 


)  LX  <  <  1 2a  t 


This  means  that  the  first  Fresnel  cone  ought  to  be  considerably  less  than  tbs 
tranirverse  dimension  of  the  perturbation.  For  a  given  wavelength  X  the  gnome  trie 

optics  approximate  formulas  hold  only  for  distances  consiierably  less  than  the 

2 

critical  length  L  B  2*1  /X.  The  absolute  magnitude  of  phase  and  amplitude 
or 

fluctuations  are,  as  one  expects,  less  affected  In  theory  by  diffraction  effects 
than  In  the  geometric  optics  approximation.  We  shall  return  to  this  matter 
after  a  statistical  treatment  of  ihs  problem. 


2.  Statistical  Treatment  of  the  Problem. 
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Up  till  nov  we  have  supposed  that  the  values  of  ths  refractive  Index 
n(x,y,s)  constitute  e  known  function  of  ooordlnatee.  In  feot,  however,  ell 
characteristics  of  the  nadlum  (including  the  refractive  index  n)  ere  constantly 
undergoing  cheotio  changes  in  tine  end  ipaoe ,  »o  that  e  proper  definition  of 
n(x,y  *,t)  can  only  be  made  on  e  statistical  besiu.  The  role  of  the  "impossible* 
value  n  *  1  will  be  that  of  the  seen  value  of  the  refractive  Index,  i.e., 

n(x,y,s)  *  1 ,  p  =  n-n. 

In  view  of  the  statistical  character  of  the  refractive  index  fluctuations 
p(x,y,s)  c  p(H)  we  lntroduoe  the  corresponding  oovariatoe  function 


BO^,  Hr,)  =  J^T^T  . 


(23) 


In  the  sequel  It  will  be  assumed  that  p(H)  is  a  stationary  process  so  that 
B(M^,  Kj)  can  be  written  as  B(r),  where  r  la  the  dlstanoe  between  the  points 
and  Kj. 

A  special  notation  for  the  real  and  imaginary  parts  of  the  kernel  of  the 


Integral  (20)  will  now  be  introduced! 


ipfa.p)-  ('Jjfu  ,p)  -  — L~rc^fL 

*r  i  _  jj a  - >  *  > '  -  rr o 

Then  the  Integral  (20)  can  be  written  s 

u  c  o)  = 


,  VI 


A 

Ab 


V  ‘  c) 


< ... )  '  I  ‘ f  ( ill  , )/.  y  C)^  ./y  < 
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fiquaring  right  and  left  aidaa  of  (24)  and  (25),  averaging  and  replacing 


the  lengthe  L,  p  and  ao  on  by  the  dlaanaionleaa  valuta  L  =  kQL,  p  *  kQp 
and  ao  on  (i.e.,  v/2*  la  unit  for  all  lengtha),  we  obtain  the  following 
•xpreaalona  for  tha  atatlatlcal  paraaatera  of  the  phaaa  and  amplitude 
fluctuation*' 


4  L  - 


SV.V1'  =  //(///  J>v4  :,f) 


26) 


1  c  a  y 
4  x  j 


L. 

I  ( 


|  l_l!j  <ft«4-4f)^i' 


where 


s  /-■  ' 


’  2 


r  -•  "■ m+'v-  f  1 

«»*  change  of  varlablea  7  =  1/2  *  ^) J  n  c  ^  -  12»  *  x  l/2  ♦  ~2)| 

^  3  ^  1  ’  (’2  “*k#*  It  poaalbl*  to  perfonr,  the  Integration  in  7  and  a.  Tor 

that  purpose  \m  mkj  eaplo7  the  faalliar  expresalon  fro*  probabl llt7  theory 


4n‘  v  V 


•  -1  <  ■*  (t ,  1) 


‘■-1 


j  !'  'P 


1..  >]%  “■  -'r , ,  , 
r.  - J'r'  = 


n  * 


»  n ' 


J  'PI 14  (<. '  f, )‘ 


.  I 


■J 


J 


("The  reeult  of  cocvolutlng  two  QauaaUn  lawa  la  again  a  Gauaaian  lav.”) 
Suiiatitutiog  flrat  In  thla  expreaalon  7^  ~  7+~  tl/^U  *i  =  *2  =  l/2  S  I 
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^  *  lOji  S2  =  -  lOg;  and  than  j1  *  J2  c(l/2^;  •1  e  *2  =  2  1  ^  *  1  *l* 

S  -  ltU;  adding  and  subtracting  the  resulting  two  aquations  and  splitting 
off  In  both  oaaaa  the  raal  part*,  on*  la  lad  to  tha  axpreaalons 


2-2  2 

whara  p  *  C,  ♦  fj  . 

It  thua  follows  from  (26)  and  (27)  that 


S'iV  0  r)'  -  j  ( f,  *-  /  ), 

(30) 

/  /  4(4  CO)V  -  1  f  J  -I  ) 

(31) 

whara 

I,- 

o  ^  „  J 

(32) 

I.  - 

[  jj  J  $(:("> >i,l,p)p(r),iric,^^j 

(33) 

r‘  =  fiLtvOV^ 

3» 
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Tbo  for.sulss  obtained  mm  to  be  rather  general  and  can  be  u»d  to  oospute 
ph*»»  and  amplitude  fluctuations  for  correlation  funotiona  B(r)  of  arbitrary 
for*,  is  an  ex**ple  lat  us  consider  tbs  following  correlation  fuactiom 

B(r)  *  n(M1)p(M2)  *  Bq  sxp 


whara  B0  *  /  and  r  is  tbs  distance  between  the  points  and  Mj. 

On,  *ay  think  that  asymptotic  laws  which  are  appropriate  to  suoh  a 
correlation  function  hold  also  (with  piecision  to  values  of  nu*srlcal  co- 
•ffloisnts)  for  any  othar  sufficiently  Mooth  oorrelatlon  function  B(r) 

which  dscreaMS  sufficiently  rapidly  at  iifinity. 

Substituting  the  function  (34)  into  fonaulas  (32)  and  (33)  enables  one 


to  explicitly  perforr.  the  integration  in  n  and^. 

Returning  now  to  the  previous  di*enslonal  quantities, 

*;  t 


find 


'V!- 


(35) 


r.  i. 


i 


( 

r 


...  i 

J 


u  /*■ 

/  'i 


(36) 


One  observes  that  the  factor  exp  r.'  inside  the  Integral 

represents  the  correlation  function  of  the  fluctuating  index  of  refreotlon 
along  the  x-axle  (i.e.f  along  the  direction  of  propagation  of  the  incident 

wave ) . 

Ft  L  >  >  the  double  integral  appearing  in  (35)  and  (36)  can  be 
•  reduced  to  a  single  integral  by  using  the  fact  that  the  function  differs 
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eeneibly  fro*  ?ero  onl  In  a  narrow  atrip  about  the  lino  ^  Tho 

correepondlng  aayaptctlc  roproaontation  la  eaaily  obtained  by  a  change  to 
tha  rariabloa  x  =  1/2 *  \2)  aad  ^  *  ^1  “  ^2*  ^  tha 

sr  -  Integration  can  be  extended  to  the  whole  real  axle  with  negligible 
error.  If  it  even  heppene  that  kQt^»  1»  then  it  la  poaalble  to  negleet 
the  luaoand  \  ^  A  with  reapeot  to  unity  I  becauae  of  the 

exponential  foot  or  the  contribution  of  tha  denominator  for 
la  aaail  In  general  J  .  Hence  there  eaaily  followa  the  following  aayaptotlc 
repreeentatlona  of  the  integrala  1^  and  1^  for  L^>>t  and  k0t^>>  1  (i.e., 

»  X/(2")  )  « 

Ij  B0  I,  *0  L  ,  (37) 


Ij  ~  d  -  B0  /,  k  n<- c  t<«" 


We  remark  that  wo  lntroduoe  the  diaenalonleaa  parameter 


(38) 


(39) 


where  It  ha  a  naturally  benn  assumed  that  L>>t  and  and  that  It  can 

in  general  be  allowed  to  take  arbitrary  values  In  computing  the  aayaptotlc 
repreeentatlona  of  1^  and  I-,. 

Froa  (30),  (31),  (37),  and  (38)  followa  the  asymptotic  expreeeiona  for 


the  mean  aquare  fluctuaticna  of  phaee  and  logarithmic  amplitude) 

7  ^  t  itK  L  (  1  p)} 


2  -  2 

We  recall  that  hare  Bq  *  n  ~  (n-n)  la  tha  Mao  square  of  tha  refractive 
index  fluctuatlona. 

Let  ua  a:w  consider  tha  limiting  casus.  for  small  valuea  of  tha  par  an  e  tar 
D  (for  D  ^  1) 

arctan  D  D,  1  -  arc  tan  D  ^  D2 


and  tha  corresponding  Kan  square  values  are 


(46) 
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°L(AJAj  ~  U  u  \  J )  (47) 

(#BI»  *1  d{a/a  )»  %  *,aa  of  the  reepeotive  quentitlea). 

Aa  on*  expect*,  theae  are  wall  known  formula a  of  gaoaatrlo  optloa  (axpraaalona 
(46)  and  (47)  a grae  with  thoaa  obtained  by  T.  A.  Araailnlkov  [8  ).  The  oaaa 

Dv<l  oorreaponda  to  comparatively  vail  value*  of  the  dlataaoa  L,  via., 
L^sLor,  where  Lcr  «=  *2/(V/2). 

For  D»1  (i.a.,  for  L>>i-Cr)  »•  h**e  arotaa  D  cc  -!■—  t  and  formula*  (40) 
and  (41)  gira 


S' 


(48) 


x  > 


it/ 


(49) 


Thua  for  large  value  a  of  L  tha  formula  for  the  phase  fluctuation! 
simplifies  to  tha  oorra a ponding  foraula  of  the  gaoaatrlc  optica  epproxima- 
tlon,  differing  froa  tha  latter  only  In  tha  numerical  factor  1/  f~2"(ln 
particular  tha  proportionality  of  «g,  to  the  equare  root  of  L  la  preserved). 
At  the  same  tine  the  dependenoe  of  tha  aaplituda  fluctuationa  00  1  change a 
aharply  upon  crone. ng  the  value  L0J.i  in  the  ease  faahion  ae  with  email 
L-valuea  the  parameter  j  increase*  proportionally  to  \?^2  aad  ia 

independent  of  wavelength  (in  agreement  with  tha  theory  of  geoaetrlo  optica). 
For  larger  raluea  of  L  the  fluctuations  of  the  logarithmic  aaplituda  follow 
the  saoje  law  aa  the  phaae  Tuctuatlonai  aiB(j^/g  )  proportional  to  Y~lT 
and  in  addition  beglna  to  dapand  on  wavelength. 


Wa  note  moreover  that,  acoordlng  to  the  above  formulae,  ia  every  oaaa 
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#ln ( A/A0)  <  V 

__  / /  ~  r  r°"  j. 

Z'  J)  +  i>  r  (  J) 

The  dependence  on  wavelength  [or  frequency)  given  by  the  theory  oan  be 

used  to  explain  the  phenomenon  of  color  changes  in  twinkling  stars. 

In  conclualon  we  shall  estimate  tba  order  of  magnitude  of  L  for 

cr 

twinkling  of  stars.  Assume  that  tba  light  wavelength  V  Is  5  x  lO-"*  an. 

and  that  4  (tba  'Internal  seal#*  of  the  atioapherlo  turbulence)  la  2  cm. 

For  theee  values  L  *  1.6  x  10^  cm.  =  1.6  kn.  77»e  value  of  the  acala  eo 
cr 

obtained  coapara a  with  tne  thickness  of  that  layer  of  the  tropoaphera  where 
one  expects  the  great  majority  of  the  tenperature  fluctuatlona  and  con¬ 
sequently  the  fluctuatlona  of  refractive  index.  Thua  oven  in  the  problem 
of  twinkling  of  atari  where  It  would  appear  that  one  ahould  find  a  very 
favorable  altuation  for  ualng  geometrlo  optloa  eethoda  (alnca  A  <<  t^)  due 
to  the  great  thlckneee  of  the  layer  one  encounters  diffraction  effecta  (in 
the  eenae  of  Fresnel). 


Qualitative  confirmation  of  tba  results  of  the  present  investigation 
can  be  gleaned  fro  a  experimental  data  on  amplitude  fluctuatlona  of  eound  and 
llgnt  wave  a  In  the  atmosphere.  Thua  naaa'uax.enta  of  the  eound  amplitude 
carried  out  by  Krasilnikov  2  ebow  that  the  fluctuation#  of  phaee  and 
logarithmic  amplitude  are  of  the  same  order  of  magnitude  but  that  the  latter 
are  so  am  what  la^s.  The  growth  of  the  logarithmic  amplitude  fluctuations 
with  increasing  frequency  ie  also  observed.  We  also  observe  in  Krasilnikov's 
results  [2  1  the  V  L- growth  of  the  logarithmic  amplitude  fluctuations  of  tbs 
eound  wove  for  larger  L  (but  not  the  law  ]?^ \  this  latter  law  of  growth 
sharply  oontraclcts  all  6t  serrations).  An  analogous  deduction  can  be  drawn 


1 
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froa  Siedentopf 's  lljaeaeureDents  of  the  intensity  fluctuations  of  light 
received  froa  distant  source i . 

Only  qualitative  deductions  froa  tbs  proposed  theory  can  b«  considered 
here.  For  a  quantitative  ooaparlson  based  on  erperiaental  results  it  is 
necessary  first  of  all  to  carry  out  calculations  for  whatever  particular 
correlation  function  is  charectarlstic  of  tbs  fluctuations  of  the  atseospherlo 
refractive  index.  The  converse  problea  le  alec  of  considerable  intereeti 
To  obteln  estimates  for  tbs  charsoterletica  of  ataoepherio  turbulenoe 
(intensity,  a  ala)  froa  statistloal  analyala  of  observed  light  intensity 
fluctuations  (twinkling  of  stars)  or  the  sound  field  characteristic  of 
•taospherlc  turbulenoe.  Consideration  of  these  problems  falls  outside  the 
soope  of  the  present  paper,  however. 


USSR  4  cede  ay  of  Saienoe 
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